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Abstract 

We show that the single-particle polaron Green's function describes SIS tun- 
nelling in cuprates, including the absence of Ohm's law at high voltages, the 
dip/hump features in the first derivative of the current, a substantial inco- 
herent spectral weight beyond quasiparticle peaks and unusual shape of the 
peaks. The theory allows us to determine the characteristic phonon frequen- 
cies, normal and superconducting gaps, impurity scattering rate, and the 

electron-phonon coupling from the tunnelling data. 
PACS numbers:74.20.Mn,74.20.-z,74.25.Jb 
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There is a number of anomalous properties of high-Tc cuprates, which seem to be intrinsic 
and universal. The isotope effect on the carrier mass unusually high values of the 
static dielectric constants [0], and the non- Fermi liquid and non-BCS spectral densities 
observed in the angle-resolved photoemission (ARPES) superconductor- insulator- 

normal metal (SIN) and superconductor- insulator-superconductor(SIS) tunnelling spectra 
|l6HT2i are among those. The isotope effect and the high values of the static dielectric 



constants as well as the optical spectroscopy |13[ suggest that electron-phonon interaction 
is more than sufficient to bind carriers and phonons into small polarons and bipolarons in 
highly polarisable cuprates [p|,|l4| . The bipolaron theory [|15| provides a natural explanation 
of the normal state pseudogap (as a half of the bipolaron binding energy, Ap) |1T6|JT7[| , NMR 
linewidth [|1^, normal state kinetics |15|, SIN tunnelling [|18| and ARPES |T9|j20|| . The theory 



accounts for two distinct energy gaps (coherent and incoherent) as observed in the recent 
tunnelling [p!0|-[T^ and Andreev reflection experiments p2|J10[| , for the parameter-free fit of 



the critical temperature , upper critical field , and specific heat [p^]. 

In this letter we present the polaron theory of SIS tunnelling in cuprates. 

Within the standard approximation the tunnelling current, I{V), between two parts 
of a superconductor separated by an insulating barrier is proportional to a convolution of 
the Fourier component of the single-hole Green's function (GF),G'(k, cu), with itself as 

/OD 
dtuQG{k, cj)53G'(p, e\V\ - uj), (1) 

K,p ~ 

where V is the voltage across the junction. 

A problem of a hole on a lattice coupled with the bosonic field of lattice vibrations has a 
solution in terms of the coherent (Glauber) states in the strong-coupling limit, A > 1, where 
the Migdal-Eliashberg theory p5|,|26| cannot be applied due to the broken translational sym- 



metry 1^^. For any type of electron-phonon interaction conserving the on-site occupation 



numbers of fermions the 1/A perturbation technique yields (at T = 0) p8|,|29 



G(k u) = ZY^ V nr=i l7(qr)P 

^' ^ i^q,^q,(2iV)'/!(^-EUi^(q.)-e(k + EUiq.)+^5)' 



where Z = exp[— (1/2A^) X]q |7(q)P]) is the matrix element of the interaction with 

phonons of the frequency uj{q) and 5— > +0 {h = c = 1). The hole energy spectrum, 
e(k) is renormalised due to familiar polaronic narrowing of the band, and (in the super- 
conducting state) also due to the interaction with the Bose-Einstein condensation (BEC) of 
bipolarons |]2T[] as 

eik) = [^kf + Af^\ (3) 

Here ^(k) = Z'E{\<.) — /i is the renormalised polaron band dispersion with the chemical 
potential fi, E{k) = X^m ■ m) is the bare (LDA) disperison in a rigid lattice. 

The mass-renormalisation exponent is 

^, ^ Em^(m)e~^'('"^ exp(-ik ■ m) 
Em^lm) exp(-ik ■ m) 

with g'^{m) = |7(q)P[l — cos(q ■ a)]. GF, Eq.(2), is exact in the extreme strong-coupling 
regime, A ^ oo (see, for example, Ref. page 294). 

Quite different from the BCS superconductor the chemical potential fi is negative in the 
bipolaronic system, so that the edge of the single-hole band is found above the chemical 
potential at = Ap. Near the edge the parabolic one-dimensional approximation for the 
oxygen hole is applied |]18|, compatible with the ARPES data [l^ 



The 'global' gap A = (Ap + A^)^/^ comprises the incoherent temperature-independent (nor- 
mal state) gap Ap, and the coherent (superconducting) gap Ac(T). The coherent gap dis- 
appears at Tc together with the condensation . 



The continuous-time Quantum-Monte-Carlo simulations ||3T| show that the 1/A pertur- 
bation result, Eq.(2), is practically exact in a wide range of the Frohlich electron-phonon 
interaction with high-frequency phonons, including the weak-coupling regime. Different 
from the canonical Migdal-Eliashberg GF there is no damping ('defasing') of low-energy 
polaronic excitations in Eq.(2) due to the electron-phonon coupling alone (because of the 



energy conservation ||3^). This coupling leads to the coherent dressing of low-energy carriers 
by phonons, which appears in GF as phonon sided bands with / > 1. On the other hand, 
the elastic scattering by impurities yields a finite life-time of the Bloch polaronic states. 
For the sake of analytical transparency we model this scattering as a constant imaginary 
self-energy, replacing i5 in Eq.(2) by a finite iV /2. In fact, the 'elastic' self-energy has been 



found explicitely as a function of energy and momentum Its energy/momentum 

dependence is essential in the subgap region of tunnelling, where it determines the value of 
the zero-bias conductance. However, it hardly plays any role in the peak region and higher 
voltages, which are of our prime interest here. 

Substiting Eq.(2) into the current, Eq.(l), and performing the intergration with respect to 
frequency and both momenta (using Eq.(5)), we obtain the tunnelling conductance, cr(V^) = 
dl/dV, 



oo 



e|V| -2A-^ 



■r=l 



E 

r'=l 



(6) 



where L [x, F] = F/(x^ + F^). To perform the remaining integrations and summations, 
we apply a model analog of the Eliashberg spectral function a'^F{uj) by replacing q-sums, 
X^q |7(q)P^(t^(q))/2A^, in Eq.(6) for the integrals {g"^ /n) ^ dujL[uj — uq,5uj]A{uj) for any 
arbitrarry function of the phonon frequency A{uj{c\)). In this way we introduce the charac- 
teristic frequency uoq of phonons strongly coupled with holes, their avarage number in the 
polaronic cloud, and their dispersion 5uj. As long as 5uj is less than cuq, we can extend the 
integration over phonon frequencies from — oo to oo and obtain 



oo „2(l+l') 

^(^) « E -jiir^ [^1^1 - 2A - (/ + /')^o, r + 5uj{i + /')] • 



«,/'=0 



(7) 



By replacing the Lorentzian in Eq.(7) with the Fourier integral, we perform the summation 
over / and /' with the final result for the conductance as 



oc / rftexp 2/e"'^'^* cos(cuot) - Ft 

JO L 



COS 



2g'^e-^^^ sin(cuot) - {e\V\ - 2A)t 



(8) 



From the isotope effect on the carrier mass, phonon densities of states, experimental 
values of the normal state pseudogap, and the residual resistivity one estimates the coupling 
strength g"^ to be of the order of 1 [|1|, the characteristic phonon frequency between 20 and 
80 meV, the phonon frequency dispersion about a few tens meV, the gap A about 30 meV, 
and the impurity scattering rate of the order of 10 meV. 

SIS conductance, Eq.(8), calculated with the parameters in this range is shown in Fig 
l,a-d for four different values of the coupling. The conductance shape is remarkably different 
from the case of BCS density of states, both s-wave and d-wave. There is no Ohm's law 
in the normal region, e\V\ > 2 A, the dip/hump features (due to phonon sided bands) are 
clearly seen already in the first derivative of the current, there is a substantial incoherent 
spectral weight beyond the quasiparticle peak for the strong coupling, > 1, and there 
is unusual shape of the quasiparticle peaks. All these features as well as the temperature 
dependence of the gap are beyond the BSC theory no matter what the symmetry of the gap 
is. However, they perfectly agree with the experimental SIS tunnelling spectra in cuprates 
||^-|TD|,|T2[ . In particular, the theory, Eq.(8) quantitatively describes one of the best tunnelling 
spetra obtained on Bi2Sr2CaCu208+5 single crystals by the break-junction technique 0, 
Fig. 2. Some excess zero-bias conductance compared with the experiment is due to our 
approximation of the elastic self-energy. The exact (energy dependent) self-energy provides 
an excellent agreement in this sub-gap region, as has been shown in Ref . . A more recent 



dynamic conductance of Bi-2212 mesas (as shown in Fig. 2 of Ref. [0]) is almost identical to 
our Fig. lb as well. The unusual shape of the main peaks (Fig. la, b) is a simple consequence 
of the (quasi) one-dimensional hole density of states near the edge of the oxygen band, 
Eq.(5). The coherent (/ = /' = 0) contribution to the current with no elastic scattering 
(r = 0) is given by 

f°° de f°° de' , ,,,,, 

SO that the conductance is a 5 function 

ao{V) (x6{e\V\-2A). (10) 



Hence, the width of the main peaks in the SIS tunnelhng, Fig. 1,2 measures directly the 
elastic scattering rate. 

The disappearance of the quasiparticle sharp peaks above Tc in Bi-cuprates has also 



been explained in the framework of the bipolaron theory [0,^. Below Tc bipolaronic Bose- 
Einstein condensation (BEC) provides an effective screening of the long-range (Coulomb) 
potential of impurities, while above Tc the scattering rate might increase by many times 



pO| . This sudden increase of F in the normal state washes out the sharp peaks from the 



tunnelling and ARPES spectra. 

Finally we would like to comment on a possible role of spin fluctuations in the tunnelling 
spectra. If they play a role, the peak-dip separation observed in ARPES and tunnelling 
should be equal to the resonance peak energy, Er, observed in the spin-flip neutron scat- 
tering . However, as discussed recently in the comprehensive review of the experimental 



constraints on the physics of cuprates [0], the peak-dip separation is nearly independent 



of Tc or the doping level, while is approximately proportional to Tc. This controversy 
as well as the direct comparison of the electron-phonon and spin fluctuation interactions 
suggest that the dip/hump features in ARPES and tunnelling arises from strong electron- 
phonon coupling (as originally proposed by one of us [|TB|). Recently, using high resolution 



ARPES data in conjuction with that from neutron, optics and local structural probes, Shen 
et al provided direct evidence for strong electron-phonon coupling being important for 
pairing. 

In conclusion, we have derived the tunnelling SIS conductance for the strongly coupled 
electrons and phonons in the (bi)polaronic regime. The theory describes SIS tunnelling 
in cuprates, including the spectral shape in the gap region, and the dip-hump incoherent 
features at higher voltages. It allows us to determine the characteristic phonon frequencies, 
the gap, impurity scattering rate, and the electron-phonon coupling from the tunnelling 
data. 
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Figure Captures 

Fig.l. SIS tunnelling conductance in the bipolaronic superconductor for different values 
of the electron-phonon coupling, g"^, and A =29 meV, ujq = 55 meV, =20 meV, F = 8.5 
meV. 

Fig. 2. Theoretical conductance of Fig. lb (solid line) compared with the tunnelling 
spetrum obtained on Bi2Sr2CaCu208+5 single crystals by the break-junction technique 

i(dots). 
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